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Using soft-collinear effective theory, the leading-log radiative electroweak corrections are written 
in a closed and analytical form for the hadronic cross section of Higgs production through vector 
boson fusion, qq —¥ qqH , one of the most promising channels for studying the Higgs boson at the 
LHC. The simple leading-log resummation is compared with a full next-to-leading-log calculation, 
and its accuracy is found to be of order 1% up to 10 TeV, i.e. comparable with the accuracy of 
PDFs. Corrections are found to be larger than predicted by one-loop fixed order approximations at 
LHC energies. The method provides a simple way of incorporating the electroweak corrections in 
software packages, improving the accuracy of simulations. 
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I. INTRODUCTION 

The recent discovery at the Large Hadron Collider 
(LHC) [H,[3 of a new boson will lead in the next years to a 
detailed study of its properties, and will require a careful 
comparison of the experimental data with the results of 
precise calculations for the Higgs sector of the Standard 
Model (SM). Vector boson fusion (VBF), qq -> qqH, is 
the second largest production channel for the Higgs bo- 
son. It is a pure electroweak process and its study is very 
important for constraining the Higgs couplings and thus 
identifying the nature of the Higgs sector [3, However, 
at the 10 Tev energy scale of LHC, even the electroweak 
radiative corrections become important 0, @] and should 
be included in computer simulations together with QCD 
corrections 

[7 14]. With a partonic center-of-mass en- 
ergy -y/i of several TeV - more than an order of magni- 
tude larger than the masses Mw,z of the gauge bosons - 
the radiative corrections contain large Sudakov logarith- 
mic terms a n L m where a = aj.,2 are the weak coupling 
constants and the logarithms L = log(s/M^ z ) emerge 
from the two different energy scales. These terms domi- 
nate the perturbative expansion and may even require re- 
summation when the fixed order perturbation expansion 
breaks down. However, for VBF the scattering ampli- 
tude is explicitly proportional to the vacuum expectation 
value (VEV), so the effective field theory operator is not 
a gauge singlet and standard resummation methods do 
not apply. That makes VBF a special interesting process 
to deal with. 

The Sudakov logarithms can be regarded as infrared 
logaritms since they diverge as Mw,z — > 0, and by using 
an effective theory they can be converted to ultravio- 
let logarithms and summed by standard renormalization 
group techniques. Quite recently the soft-collinear effec- 
tive theory (SCET) [l5l - [l~8| has been shown to provide a 
simple way to obtain the sum of the series of leading-logs 
(LL) a n L n+1 , next-to-leading-logs (NLL) a n U\ next-to- 
next-to-leading-logs (NNLL) a n £ n -\ etc. pjl^. In 
the effective theory the single terms contributing to the 



scattering amplitude are multiplied by the general fac- 
tor El 



U(lH,(i h ) = cxp 



x exp 



x expF(a(fj, h )) 



D (a(w))+-Di(a(w))log^ 

AK/x)) log ^- +£(«(//)) 



(1) 



where Hi ~ Mz is the low energy scale and [ih ~ \fs is 
the high energy scale. The coefficient F is the high scale 
matching coefficient, D and Di are the low scale match- 
ing coefficients, B is the non-cusp anomalous dimension 
and A is the coefficient of the cusp anomalous dimension. 
The LL series is summed by the one-loop cusp anoma- 
lous dimension; the sum of the NLL series requires the 
two-loop cusp anomalous dimension, the one-loop non- 
cusp anomalous dimension, and the one-loop low scale 
matching coefficient D\ ; the NNLL series is given by the 
the three-loop cusp, two-loop B and D\, one- loop D 
and F, etc. In fact the exponentiated form of Eq.(JTJ) 
only requires the inclusion of electroweak corrections at 
low orders, while the unexponentiated form of fixed-order 
calculations would require electroweak corrections of any 
order for achieving the same accuracy. In some recent 
papers 0, [24[ the method as been used for calculat- 
ing the electroweak corrections to Higgs production via 
VBF. Numerical results at NLL order were obtained for 
the cross section, inclu ding the effect of parton distribu- 
tion functions (PDFs) [2J| . Most of the corrections were 
obtained in analytical form by SCET and might be easily 
included in the other software packages that have been 
developed, without the need of tedious loop calculations. 

While PDFs imply a 1% error, many of the retained 
terms at NLL order are very small, less than 1%, and 
could be neglected in order to speed up the numerical 
integration of the cross section. In this paper the simpler 
LL order resummation is compared with the full NLL cal- 
culation for Higgs production via VBF, and the accuracy 
of the two approximations is discussed. The LL calcu- 
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lation of the cross section is found to deviate from the 
NLL result by less than 1% up to 10 TeV center-of-mass 
energy. While the present study gives a useful evalua- 
tion of the accuracy in the effective theory resummation, 
it provides a very simple and fast analytical way for in- 
cluding the electroweak radiative corrections in software 
packages. 

In fact at NLL order the one-loop anomalous dimen- 
sion turns out to be a 10 x 10 matrix in the operator 
basis [Hj], and the running integral in Eq.fll} requires a 
numerical computation. Moreover two-loop beta func- 
tions are required for the running of the couplings at the 
same order, yielding coupled equations that again must 
be solved by a numerical routine. On the other hand, 
at LL order the only one-loop term that is required is 
the cusp anomalous dimension, and by use of the simple 
uncoupled one-loop beta functions the running integral 
in Eq.(TT]) is analytical, yelding a diagonal correction fac- 
tor for the differential cross section. Moreover, an exact 
cancellation of terms yields the same correction factor 
that would be found for the quark scattering qq — > qq, 
plus Higgs rescattering and wave function renormaliza- 
tion terms. NLL terms are still shown to be relevant for 
a full description of the dependence on scattering angles 
of the differential cross section, but such small depen- 
dence is averaged in the integrated cross section. 

The paper is organized as follows: in Sec. [II] the kine- 
matics of the VBF process is described and some de- 
tails on the integration of the cross section are reported; 
the general calculation framework is described in Sec. IIIII 
where the operator basis set is defined for the effective 
theory and matched onto the low scale and high scale 
gauge theories; in Sec. |TV]tlie running of the Wilson co- 
efficients is discussed and the anomalous dimension is 
shown to take a simple diagonal form at LL order; ex- 
plicit analytical expressions are derived for the running 
integral in Sec. [V] and the numerical results at LL or- 
der and NLL order ar compared with other fixed order 
perturbative calculations in Sec. EH 



II. KINEMATICS OF VBF AND CROSS 
SECTION 

The tree-level diagrams for the VBF process qq —> qqH 
are shown in Fig.l. The two outgoing quarks produce for- 
ward jets with a large rapidity gap, which characterizes 




FIG. 1: The tree-level t, u and s-channel diagrams for elec- 
troweak Higgs production. 



the VBF channel. This allows one to suppress reducible 
and irreducible (i.e. gg — > H + 2 jets) backgrounds us- 
ing selection criteria called VBF cuts. These consist of 
tagging the two forward jets and sometimes a veto on cen- 
tral jets, see e.g. Refs. [8|, l25H27j . The momenta of the 
incoming quarks (beams) are denoted by p\ and pi and 
the momenta of the outgoing quarks (jets) are denoted 
by j»3 and p^. The Higgs boson momentum is ph and we 
will assume that it is produced on-shell, pi = M%. We 
will neglect quark masses, p\ = p i, = p\ = p\ = 0, and 
follow the conventions of Ref. |23[ by taking all momenta 
to be incoming. The kinematic configuration is fully de- 
termined by the following generalized Mandelstam vari- 
ables: 

S = {pi + p 2 ) 2 , t 3 = {pi + p 3 ) 2 , t 4 = (p 2 + Pif , 
U3 = {P2+P3) 2 , U4 = (Pi +P4) 2 ■ (2) 

In SCET, each collinear direction has a set of two light- 
cone reference vectors rii = ±(l,n.;), hi = ±(1,— n«) 
associated with it, where in our conventions we take 
the plus (minus) sign for incoming (outgoing) parti- 
cles. For the quarks we can simply take rii = ±pi/p®. 
For the Higg s we have t o take its mass into account, 
n„ = - Vh /^{plf~Ml 

At hadron colliders, such as the Tevatron and LHC, 
the colliding quarks (i = 1, 2) move along the beam axis 
and carry a fraction Xi of the hadron momenta, 

p»= Xl E cm ^-, <= (1,0,0, ±1), (3) 

where E cm is the center-of-mass energy. In the cross sec- 
tion, the momentum fractions Xi are integrated over. The 
distribution of momentum fractions carried by the quarks 
in the proton is described by PDFs. We parametrize the 
momenta of the outgoing quarks by 

P3 = --E 3 (l,sin03,O,cos0 3 ) , 

Pi = —Eiil, sin #4 cos if, sin #4 sin (p, cos #4) , (4) 

where Ei > is energy and 0i the angle with the beam 
axis of particle i. The 9i are related to (pseudo)rapidities 
r\i by tan(#i/2) = exp(— rji). The azimuthal angle be- 
tween the particles is given by ip. Momentum conser- 
vation fixes the momentum of the Higgs and yields one 
additional constraint. 

In an effective theory, to obtain the electroweak cross 
section, we need to square the matrix elements of the ef- 
fective operators, sum over quark helicities, flavors, chan- 
nels and integrate over phase space. Schematically, 

1 f dxi dx 2 1 \ - \ - \ - 

cm J 1 z hoi. Hav. chan. 

x ]T Mx 1 )f,(x 2 ) C* x C Y (d x )*(6y) , (5) 

where Ox are the effective operators and fi(x) denotes 
the PDF of flavor i evaluated at momentum fraction x. 
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An identical particle factor of 1/2 (not shown) must be 
included for symmetric phase-space integrations. The 
Wilson coefficients C are obtained in terms of tree-level 
high scale coefficients C by matching at the high scale, 
running down to the low scale and matching at the low 
scale according to Eq. (P) , 



C(ni) = U(huHh) ■ C(hh) 



(6) 



which must be regarded as a matrix product in the op- 
erator basis set. The three-body phase space d$3 reads 



/ d$ 3 = / n 

J J i=3,4 



d 3 pi d^ph 
(2tt) 3 2E 1 (2tt) 3 



e(- P ° h )6(pl-M* 



x (2tt) 4 5 4 ( Pi + p 2 + p 3 + Pi + Ph) 
1 



2 6 tH 



[EidEid(coBOi)]dtp0 



i=3,4 



x (5 



4 

£ 

3=1 



ft 



— Mi 



2 6 7T 4 



d(cos 3 ) d(cos 84) dip dE- c 



(7) 



In the first step we carry out the pu integral using the mo- 
mentum conserving delta function, and write pi in terms 
of the spherical coordinates in Eq. (QJ. In the second step 
we use the on-shell condition for the Higgs to perform the 
£4 integral, yielding £4 = /1//2 with 



/1 = xix 2 E cm - M h - [xi + x 2 + (x 2 - xi) cos 9 3 ]E cin E 3 

f% = [%l + x 2 + (x 2 — X\) COS #4] E'en — 2[1 — COS 63 COS 04 

— cos ^ sin 03 sin04]i?3 . (8) 

The remaining integrals in Eq. must be carried out 
numerically. Their boundary conditions are 



< 03,04 < 7T, 

Mi 

< a?i < 1 , 



B 2 

< £3 < 



< < 2vr , 

< x 2 < 1 , 



xix 2 E* 



Ml 



[Xl + X 2 + (X2 — Xi) COS0 3 ]i? c 



(9) 



where we included the bounds on momentum fractions of 
the PDFs. The boundaries will be further restricted by 
any cuts we impose on the final state. 



III. EFFECTIVE THEORY AND OPERATOR 
BASIS SET 

Electroweak corrections can be obtained by SCET in 
the framework of Ref. [19i. l22|. T he extension to the VBF 



process was derived in Ref. 1231. |24|] . and the explicit ex- 
pressions are reported in Ref. 24[ . Here we only consider 
the SM gauge group. Extensions like the minimal left- 
right symmetric gauge group [28l43ll | will be the subject 
of an other paper. 



The first step consists of matching onto SCET at a high 
scale fih ~ y/s. This matching can be done in the un- 
broken gauge theory, since symmetry breaking effects are 
suppressed. Using the renormalization group evolution, 
one then runs the effective theory operators down to a 
low scale \n ~ Mz- At the low-scale, the W and Z boson 
are integrated out and one matches onto a SU(3) x U(l) 
effective theory, which only contains gluons and photons. 
The effects of SU(2) x U(l) symmetry breaking only en- 
ter in this low-scale matching. 

At the high scale, we can ignore the effect of elec- 
troweak symmetry breaking and work in the unbroken 
phase of SU(2) x U(l), since the partonic center-of- 
mass energy y/s is large compared to the VEV. The ba- 
sis of SCET operators for vector boson fusion is given 

by mm 



1A ,B = 1 A ,0 1 B , 
2A , B ,c = O a 2 A ,O a 2 O a B ,O a 2 O a c , 

3A .B = 3 A , 3 O b , 
G 4A ,B,c = OtO A ,OlO B ,OtO%. 



(10) 



The Higgs sector is described by 0\ , . . . , O4 
Oi = $j> + 4>l&h : 



O a 2 = & h T a <t> 



drna 



T a <5> h 



3 = * 



(11) 



where $/, = Wl. 



denotes the collinear scalar doublet 



<f) nh that will produce the Higgs, plus the corresponding 
collinear Wilson line W nh . The field <po denotes a soft 
scalar that will attain a VEV in the broken phase. The 
basis of operators for the quarks is 

A = *37 M r a * 1 *47M :ra *2 , 

O b = ^#37^*1*47^*2 , 

0\ = *37 M r a # 1 *47M*2 , 
0% = §37^*1 *47 M T °*2 , 



O a c = ie abc #37 M T fc #i*47 A1 T c # 2 , 



(12) 



where the subscript i = 1 , . . . , 4 on the field labels the 
particle with momentum pi . All these fermion fields con- 
tain collinear Wilson lines ^ = W£.i/) ni , as required by 
collinear gauge invariance. In Eq. (1121) . is a fermion 
doublet (singlet) if it is left-handed (right-handed). We 
suppressed the projectors Pr.l = (1 ± 7s)/2 to keep our 
notation general and allow for both left- and right-handed 
quarks. We can consider each helicity separately and 
combine the contributions at the end. Operators O a and 
Oq are only well-defined if all quarks are left-handed, and 
O a A and Og are only allowed if at least one of the quarks 
is left-handed. 

The basis of operators in Eq. (TT2]) is not complete, be- 
cause it assumes that the incoming particles are quarks 
and because it only suffices for the t-channel contribution. 
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If one (or both) of the incoming particles is an anti-quark, 
we can simply obtain the corresponding basis and expres- 
sions by interchanging 1 <H> 3 or (and) 2 O 4. Similarly, 
to switch to the it-channel we interchange 3f*4. We can 
thus obtain the results for the other channels by simply 
making replacements and will therefore restrict the dis- 
cussion to the t-channel with incoming quarks. 

At LL and NLL order the tree-level high-scale match- 
ing suffices, and can be done in the unbroken phase of 
the electroweak gauge theory. The relevant terms in the 
Lagrangian that couple the scalar doublet to the gauge 
fields are 



+ g l g 2 B tl A a u &T a <5>. 



(13) 



Finally, let us consider the low energy matching. The 
tree-level result suffices at LL order. The one-loop calcu- 
lation was required at NLL order and was also reported 
in Ref. (24|. At low energies the effects of electroweak 
symmetry breaking need to be taken into account. We 
match onto a basis of operators in the broken phase of 
the gauge group. 

We first perform the tree-level matching for the Higgs 
part of the operator, given in Eq. (fTTj) . At tree level, the 
soft scalar field simply attains a VEV and the collinear 
scalar field produces a Higgs, leading to 

{O x , O a 2 , 3 , 0%) -► {1, 0, 0, ~S a3 /2}vh (17) 

where v is the VEV and h is the Higgs field. 

The quark part in Eq. (fT2|) gets matched onto the set 



The full-theory diagrams were shown in Fig.l. Depend- 
ing on the parton types, only some of the diagrams con- 
tribute. Phenomenologically, vector boson fusion is most 
interesting when the jets are in the forward direction, 
for which the s-channel contribution is suppressed. We 
therefore do not include the s-channel. Matching the t- 
channel diagram for incoming quarks onto the operators 
in Eq. (fTUl), yields 



Cia{h) = . . vm , Cibw = -, , Y1Y2 ■ 



2t 3 t4 

c 4A w = — — mYi 

t 3 t 4 

- r-7 ¥ lV2 



(14) 



The other Wilson coefficients vanish at tree level. Here, 
the variable rji, is defined as rji = 1 if the particle with 
label i is left-handed and rji = if the particle is right- 
handed. The hypercharge Y — 1/6 for left-doublets, Y — 
2/3 for right-handed up- type quarks and Y = —1/3 for 
right-handed down-type quarks. 

The running of the Wilson coefficients is described by 
the RG equation 



d/z 



(15) 



where 7 is the anomalous dimension for the operators. 
It should be noted that this is a matrix equation, cor- 
responding to the 10 operators in Eq. (fT0|) . Using the 
notation of Eq.([T]) the anomalous dimension can be writ- 
ten as the sum of cusp and non-cusp terms 



O a 




O b 


= u^uid-il^ 


Oc 


= d 3 7 Prf l«47M u 2 


o D 


= d^dxd^^ 




= d^uiu^^ 


6 F 


= -0 3 7 M did47/i u 2 



(18) 



Here u and d denote up and down-type fields. For op- 
erators Oa , ■ ■ ■ , Od , each pair of fields ^j^ip is a pair of 
left-handed or right-handed fields, whereas in Oe and Of 
all fields are left handed. We match {Oia,b, 04A,b,c} 
onto the quark operators {Oa 7 ■ ■ ■ , Of} in Eq. (TTSj) . ig- 
noring O2 and O3 which vanish at tree level. This is 
described by a 6 x 5 matching matrix , 



Oi 



such that 



vh^0 3 Rf : 



(19) 



(20) 



The matrix is given by 





/ 1 


3 


-1 


-1 






-1 


3 


-1 


1 





1 


-1 


3 


1 


-1 





4 


1 


3 


1 


1 







2 











-1 




V2 











1/ 



(21) 



A(a(jO)]og^ 



B(a(n)) , 



(16) 



We only need the one-loop cusp anomalous dimension 
at LL order, while two-loop cusp and one-loop non-cusp 
terms were required at NLL order and were reported in 
detail in Ref. 0. 



IV. ANOMALOUS DIMENSION AND 
RUNNING 

At LL order, one-loop corrections enter through the 
cusp anomalous dimension. For VBF the one-loop terms 
contributing to the anomalous dimension were derived in 
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Ref. [23J and reported in detail in Ref. [24[ . They can be 
written as the sum of three terms 



7 = 7c + is + 1- 



(22) 



The first term 7c is the collinear anomalous dimen- 
sion, it is diagonal and contains the large logarithms 
that contribute to the cusp anomalous dimension, plus 
wave function renormalization terms. Resummation at 
LL order requires the inclusion of these one-loop terms, 
and they are easily obtained by the sum of single-particle 
collinear terms [12, 112] for all the external particles 

7c = £>c+7c (23) 

i 

As reported in Ref. [24|, with the shorthand 



Li = log 



n% - Pi 



the quark one-loop collinear terms Jq read 



lb 



3«2 
47T 



Y- 



9«2 
16tT ? 



3ai 

47T 



Y 



while the Higgs one-loop collinear term 7^ is 



7c 



3a 2 

47T 



Ol 
47T 



if, 



3«2 
47T 



a 1 

47T 



16tt 2 



(24) 



(25) 



(26) 



where we included the contribution from the top Yukawa 
yt to the Higgs wave function renormalization. 

The second term 75 in Eq. (|2"2l is the soft anomalous di- 
mension that is obtained by summing over the soft func- 
tions [22], |32| , and is reported in Ref. 24] in some detail. 
This soft term does not contain any large logarithm and 
does not contribute to the cusp anomalous dimension. It 
only depends on the scattering angles of the external par- 
ticles, and its effect largely cancels in the integration of 
the cross section. At LL order this term can be neglected. 

The third term 7 in Eq. (|22[) is a specific new contri- 
bution occurring in the VBF process [23| . and can be 
written as the sum of SU(2) and U(l) parts plus a term 
7 A arising from the rescattering of the Higgs boson. The 
first two parts contain cusp diagonal terms proportional 
to log /i, and smaller non-cusp off-diagonal terms that do 
not depend on the scale fi. At LL order we only need 
to retain the diagonal cusp terms and the rescattering 
term that is non-cusp but is large because of the large 
self-coupling A of the Higgs boson. We also include the 
diagonal wave function renormalization terms. Accord- 
ing to Ref. [lij], the cusp parts of the SU(2) contribution 
read 



-SC7(2 
hA,U 

for the coefficients Cia, C\b 



7ii,ii = -^[2L fc + l] 



~S£7(2) 



Q 2 

8^ 



[2£ h -4r?i(ii+£ 3 ) + l] 



(27) 



(28) 



for the coefficient C4A and 

i!b (2) = ~ 4r? 2 (L 2 + U) + 1] (29) 

for the coefficient Cab- The cusp part of the U(l) con- 
tribution is 



7 U(1) = -^[2^ + l], 



(30) 



and the Higgs rescattering contribution is given by 



7 



4tt 2 



(31) 



where c° = {3,3,1,1} for the coefficient 
{Cia, Cib, C4A, Cab}- 

Summing all the terms, we find that the dependence 
on L/j cancels (up to non-cusp terms like the difference 
Lh — Li that can be neglected at LL order), yielding the 
following simple result for the total anomalous dimension 



^5E^ + 7E^ + ^ ( 32 ) 



where the sum is over the four external quarks. Here 
7„ c is the sum of the retained diagonal non-cusp terms, 
including constant wave function renormalization and 
Higgs rescattering terms 



In 



9a 2 
"16tt 



3ai 

47T 



E^ 2 



16tt 2 



+7 A 



(33) 

where d° = {2,2,10/9,10/9} for the coefficient 
{C\A,CiB,CiA,CiB}- Some of these terms are small, 
but they are constant, do not depend on scattering an- 
gles, and their weight might sum up in the integration 
of the cross section. All other non-cusp terms have been 
neglected. 

The simple result of Eq.([3"2"j) says that the one- loop 
cusp anomalous dimension is the same that we would 
obtain by the sum of the collinear terms for the external 
quarks, neglecting the Higgs particle. However the Higgs 
momentum would affect the kinematic of the quarks any- 
way. Moreover the Higgs rescattering and wave function 
renormalization terms are not small and have been in- 
cluded in the non-cusp part 7„ c of the anomalous dimen- 
sion. 



V. RADIATIVE CORRECTIONS AT LL ORDER 

At LL order the running of the coupling constant 
can be evaluated by uncoupled one-loop beta functions. 
That, together with the simple diagonal form of the one- 
loop cusp anomalous dimension in Eq. (l32[) . allows for a 
fully analytical evaluation of the electroweak radiative 
corrections. 
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By insertion of Eq. (fT9l) in Eq.([T]), the low scale Wilson 
coefficients follow from Eq.© that now reads 



<7y(«) = vJ2R { yx \ U A xO*i^h) \ Cx(Mh) (34) 

where the product runs over the four couplings a — 
a>i,a2,a y ,a\, having defined a y — yf/(47r) and a\ — 
A/(47r), and the exponentiated running factors follow by 
integration of the corresponding anomalous dimension 



exp ■ 



d/i 



(35) 



where the one-loop coefficients are well known [3 5 



h = 12^' 



b 2 =- 



19 

127T ' 



bx =- 



n 
9 

47T 



(41) 



By a simple integration, the running factors in Eq. (|35[) 
take the explicit general form 



Here 7^ is the term proportional to the coupling a in 
the anomalous dimension for the Wilson coefficient Cx- 
By inspection of Eq. (l32j) , we see that such terms can be 
written as functions of t — log fi 



(a) 

Tx 



Ctct + jjpfp) 



(36) 



where for each coupling a there is a different coefficient 
a a , while JF-^ (p) is a function of the external momenta 
p = (px,P2,Ps,P4) that in general also depends on the 
chosen coefficient Cx- The explicit expressions of the 
functions are 



3 

4-7T 



m fog(^ - Pi) - 



-in 
3 

4^r ' 



16tt 



E^ 2 



.0 



^ a) (p)=-, 

7T 

while the coefficients are 

a x = - 

a.2 = - 



(37) 



IT ^ 

i 



7T 

_3_ 

47T 



a\ =0 , 



(38) 



In all these definitions the sums are over the external 
quarks. The couplings also have an implicit dependence 
on i, dictated by the one-loop beta-functions 



doj 



These can be easily integrated, yelding 

' Oi(tl) 



Cti(th) 



= 1 - biai(ti)(t h - ti) 



(39) 



(40) 



■A> x '{null*) = — 



a{ti) 
a{t h ) 



1 



b a \b a a(ti) 



+ a a t l +F i x ) (p) 



(42) 



and their product in Eq. (|34)) gives an analytical expres- 
sion for the electroweak radiative correction at LL order. 
The cross section follows by insertion of the low scale Wil- 
son coefficients from Eq. (|Mf in the phase space integral 
of Eq.© and integrating by PDFs. The closed analytical 
form of Eq. (|42j) greatly speeds up the numerical integra- 
tion of the cross section at LL order compared with NLL 
order. We discuss the accuracy of the two orders in the 
next section by a direct comparison. 



VI. LL VS. 



NLL ORDER: 
RESULTS 



NUMERICAL 



In this section the accuracy of the approximation is 
tested by a comparison between LL and NLL orders. All 
the details of the numerical integration are kept exactly 
the same as reported in Ref. [24j where the cross sec- 
tion was evaluated at NLL order. We summarize them 
briefly. The low energy matching scale is chosen to be 
lii = Mz ■ The couplings and parameters of the standard 
model have been set at the electroweak scale = Mz 
according to the data of Ref. [3J| . The Higgs mass is as- 
sumed to be Mh = 125 GeV. The high scale fih is set 
at the larger value between Mz and the geometric av- 
erage /i 2 = \ft$t±. At this scale the sum of logarithmic 
terms ln(— t^/fi 2 ) + ln(— t^/fi 2 ) reaches its minimum in 
the one-loop matching. As discussed in Ref. [24| this 
choice has the merit of stopping the running whenever 
one of the Mandelstam variables is too small, while keep- 
ing the neglected one-loop matching terms as small as 
possible. On the other hand the sensitivity of the result 
to the choice of the high scale was shown to be small and 
comparable to the sensitivity of the standard tree-level 
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FIG. 2: Single terms contibuting to the tree-level cross section 
as a function of the center-of-mass energy E cul . Cuts are: 
p T > 20 GeV, 6*3 > 10°, <9 4 < 170°. 
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FIG. 3: The K-factor K = oew /&tree for the process u^d^ — > 
d,LUL as a function of 8 = #3 = 180° — 64, at tp = 90°, x\ = 
x 2 = 0.5, E crn = 14 TeV and E z = Zi^-Ecm/ 2 - The result 
at LL order (solid line) is compared with the NLL order data 
points (squares). 



cross section. 1 We use CTEQ6 PDFs|39| and neglect the 
very small contribution of t and b quarks. Moreover we 
neglect the s-channel contribution and interference terms 
that are known to be small with VBF cuts. The masses 
of the vector gauge bosons are restored in the denomi- 
nators of the coefficients in Eq. (fT4")) : while that is not 
so relevant when relative cross sections are reported with 
large cuts, the effect of masses becomes crucial when the 
total cross section is evaluated by integrating over the 
whole phase space. Hereafter, in order to compare with 
the NLL calculation of Ref . 24J , we adopt the same cuts 
on angles and transverse moment: 6*3 > 10°, 64 < 170°, 
Pt > 20 GeV. No QCD corrections have been included in 
both calculations, and only the virtual electroweak cor- 
rections are considered by the method. 

The single terms contributing to the cross section at 
tree-level are reported in Fig. 2. The ud — ¥ du process 
is the dominating one, followed by the other Left-Left 
terms. The Left-Right terms are two order of magnitude 
smaller, while the Right-Right terms are negligible and 
have not been included in the calculation. A sum over 
generations of the up- and down-type quarks is included 
through the PDFs in the terms of Fig. 2. 

As a first comparison we study the single terms in the 
phase space, and for each of them we calculate the K- 
factor which is defined as the ratio between the cross 
section with electroweak corrections on, and the tree-level 



At variance with the variational approaches to the electroweak 
sector [35l437| . the principle of minimum sensitivity here does not 
apply even if /i h is not a physical scale. In fact the sensitivity 
only measures the dependence on fi^ of the neglected terms, and 
this dependence might be high even when the neglected terms 
are small. The best choice of fi^ should make the omitted terms 
small, as for the method of minimal variance |38H. 
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FIG. 4: The same as Fig.3 for the process Uhdh — > Uhdh- 
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FIG. 5: The same as Fig.3 for the process di,d_L — ► dLdL- 
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FIG. 6: The same as Fig.3 for the process ulUl — ► ulUl- 
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FIG. 8: The same as Fig.3 for the process UlUr — >■ ulUr. 
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FIG. 7: The same as Fig.3 for the process urcLl urcLl- 
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FIG. 9: The same as Fig.3 for the process cLlcLr — > (IlcIr. 



cross section without any radiative correction 



K = 



(43) 



where the differential cross sections are evaluated for 
fixed values of x\,X2, E3 and E cm , and for a given set 
of angles. Here the simple LL calculation is compared 
with the NLL result of Ref. 0. In Figs. 3-10 the PC- 
factor is reported as a function of 9 = 9 3 with 9 4 = tt — 9% 
and ip — 7r/2 at a typical set of parameters: E cm = 14 
TeV, x\ — X2 — 0.5 and E% = x\X2E cla /2. The agree- 
ment of LL and NLL results is very good at small angles. 
As expected, at LL order the dependence on angles is 
reduced in comparison with the NLL result, because the 
neglected terms have a larger dependence on the scatter- 
ing angles. However for 9 < 90°, in the physically rele- 
vant phase-space region, the difference is small and the 
LL result (solid line) interpolate the NLL data (squares) 
quite well. The dependence on the azimuthal angle <f> is 
very small and negligible for 9 < 90°, in perfect agree- 
ment with the NLL order. 

We find a better agreement for the Left- Right processes 
(Fig. 7-10) where the radiative corrections are larger, 



co 0.6 
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FIG. 10: The same as Fig.3 for the process ul<1r — > ul<1r. 



than for Left-Left processes (Fig. 3-6) where the cor- 
rections are smaller. Thus the very large suppression of 
the Left-Right processes is mainly due to the role played 
by the collinear anomalous dimension. It is remarkable 
that such large suppression has no relevant effects on the 
total cross section which is dominated by the Left-Left 
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tree as a function of the center-of-mass en- 



FIG. 11: 

cttree]/ '<Jtr 

ergy E cm evaluated by Eq. ([5} at LL order (filled squares), 
compared with the NLL order of Ref. [24[ (filled circles) and 
with the output of HAWK (open circles) . Cuts are: pr > 20 
GeV, 6» 3 > 10°, 8 4 < 170°. 



processes. 

The single terms are integrated by PDFs over the phase 
space and summed up, yielding the cross section (Jew- 
By comparison with the tree level cross section, the rel- 
ative electroweak correction is defined as the ratio 



A(T _ (JEW ~ &tree 
<J (J tree 



and is displayed in Fig. 11 as a function of the center-of- 
mass energy. For comparison the NLL result of Ref. [HJ 
is also reported in Fig. 11 together with the output of the 
code HAWK [5, 6] that is based on a fixed order pertur- 
bativc calculation. 



(44) 



While differences are negligible below 2 TeV, at large 
energies LL and NLL corrections grow faster than pre- 
dicted by fixed order one-loop calculations, reaching 9% 
at the full LHC energy E cm = 14 TeV, to be compared 
with 5% predicted by HAWK. At LL order the correc- 
tion is a bit smaller with respect to NLL order, but the 
difference is less than 1% up to E cm = 10 TeV. Thus the 
simple LL calculation provides an analytical electroweak 
correction that seems to be more reliable than fixed order 
calculations at high energies. While fixed-order pertur- 
bative calculations might miss part of the correction at 
the LHC energy scale, the simple LL resummation con- 
tains the main terms, and gives integrated corrections 
that are as accurate as PDFs. Insertion of the simple an- 
alytical result of Eqs. (f34]) .(j42 j) in software packages would 
be straightforward, and would increase the accuracy of 
simulations. 
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